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Abstract 

Based on an analytical technique using a unitary transformation and the variational method, 
we study the chiral order parameter in the Schwinger model in the lattice formalism with Kogut- 
Susskind fermions. The fermion condensate (vl/ "f) for any coupling constant and fermion mass 
are calculated. Chiral symmetry is shown to be broken in the massless limit and good scaling 
behavior is obtained. 
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1 Introduction 



In the field of lattice gauge theory (LGT), there are two major directions of studies. One is to 
create new techniques such that Monte-Carlo (MC) simulation for systems with fermions so as to 
get exact results efficiently. The other is to investigate the universality of lattice action and the 
existence of the continuum limit of LGT. 

The Schwinger model|jl], which describes 1+1 dimensional QED, is superrenomalizable and 
exactly solvable[||] - Previous investigations of the model reveal some important properties of 
QCD such as quark confinement, chiral symmetry breaking and U{\) problem. LGT with fermions 
in 1+1 dimensions is easier to study than in higher dimensions. The study of lattice formalism of the 
Schwinger model is a good test of the two major directions of LGT. Before we apply the technique 
to more realistic systems, the studies of the Schwinger model can give us a better understanding 
of LGT. 

The investigations of the Schwinger model have been carried out by MC simulation and by 
other analytic approach - [0]. Although MC simulation gives several exciting results, analytic 
methods are still necessary to realize the quatitative picture of its lattice formalism. However, even 
the finite-lattice technique |T5| - |l^, one of the most successful analytic methods, is restricted to 
the strong coupling regime. New analytic approaches still need to be explored. 

Recently, a technique with naive fermion formalism||l^ is developed. It contains a unitary 
transformation and the variational method such that the Hamiltonian can be diagonalized and 
reliable physical informations can be extracted from the Hamiltonian LGT. Our study is concen- 
trated on the Kogut-Susskind formalism of this approach for the Schwinger model. The chiral order 
parameter {^i ^) is calculated and its scaling behavior is investigated. 

This paper is organized as follows. In Section 2, we give a brief review of the Kogut-Susskind 
formalism for the Schwinger model. The unitary transformation and the variational method are 
introduced in Section 3. Finally, the results and discussion are given in Section 4. 

2 The Schwinger model on a lattice 

We consider the Schwinger model on a lattice and use the Hamiltonian formalism with Kogut- 
Susskind fermions. In this formalism, the upper and lower components of the fermion spinor are 
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put on different lattice sites. We shall work in tfie temporal gauge and the Hamiltonian reads 



^ = -f E [x'{x)U{x, X + l)x(x + 1) + h.c] + mo E (-l)"x+(^)x(x) + |^ ^ L\x), (1) 

X XX 

where x"'"(x) is a fermionic raising (lowering) operator at site x, U{x,x + 1) is the gauge field 

on the link between site x and x + 1, a is the lattice spacing, mo is the bare mass of the fermion 
and g is the dimensionless bare coupling constant related to the charge hy g = ae. U{x,x + 1) can 
be written as 

U{x,x + 1) = e^""^(^) = e^^(^) = 6+(x). (2) 

The dimensionless operator L{x) which generates rotation in 9{x) is the canonical conjugate to 
e{x): 



Specifying \n) as the eigenstate of L{x), we see that h'^{x) {h{x)) raises (lowers) the boson number 
by one: 

6+ In) = In + 1) , 

(4) 

h \n) = |n — 1) . 

We also need the commutation relations between the boson operators and the fermion operators: 



[x+(x),6+(x)] = [x{x),b+{x)\ = [x+{x)Mx)\ = [x{x)Mx)]=Q- (5) 
The fermion operators % and x"^ satisfy the usual algebra: 

{x+(^), x{y)) = Sx,y, {x-^ix), x+(y)} = {xW, x(y)} = o. (6) 

If we consider the dimensionless parameter g^ to be large (strong coupling), the first term 
in the Hamiltonian can be considered as a perturbed term. The second and the last terms are 
the unperturbed terms. The ground state |0) of the unperturbed Hamiltonian has the following 
properties: 
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|0) = 0, for all x, 

0, for even x 

X+(x)x(x) |0) = { 

|0) , for odd X. 

Because odd lattice sites are occupied by an antiparticle, this state corresponds to a "filled Dirac 
sea." 

The gauge invariant states are specified by fermion occupation numbers and by boson numbers. 
They are denoted by 

m=\l{x'{x)y^^\h+{x))<^^o), (8) 

X 

where |0) is the "filled Dirac sea," and /(x), s{x) are the numbers of fermions and bosons at location 
X, respectively. 

3 Unitary transformation and the variational method 

We use the unitary transformation: 

H = exp{-ieSf)Hex-p{ieSf), (9) 
to diagonalize the Hamiltonian. The physical vacuum is assumed to be 

\n) = ex.p{ieSf) |o) , (10) 

where 



^/ = ^ E (-1)" [x^i^Pix, X + l)x{x + 1) - x^{x + l)U+{x, X + l)x{x)] , 

X 

and 9 is the variational parameter. 

The physical vacuum energy is then given by 

Eq = (0| h' |0) = (0| e^pi-ieSf)Hexp{ieSf) |0) . (11) 

Defining 



H = exp{-ieSf)Hexp{i9Sf) = Hj^ + H^ + Hg, (12) 
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where 



Hj^ = eM-iOSf) (-^ E [x'-{x)U{x,x + l)x{x + 1) + h.c.]) exp(z05/) 

H'^ = exp{-i9Sf) mo E {-irx+{x)x{x) expiiOSf) (13) 



and using the Baker-HausdorfF formula 

e^Fe-^ = F + [A,F] + ^^ [A, [A, F]] + i [A, [A, [A, F]]] + , 

and the commutation relations, we can work out H'f,, H'^, and Hg. 

We first consider H'^. In the expansion of the Baker-Hausdorff formula, we find that, because 
there is one link in 5/ and none in H'^, there appear odd numbers of U, combinations for the 
even terms. 

(0| C/"^(C/+)" |0) = if m + n is odd. (14) 

We therefore conclude that 

(0| H'^ |0) = J^.i-V2er''rno{-irO E (-1)^ (0| X+(x)x(x) |0) 

n=0 ^ ' 



= -mo^Jo{2V2d), 

where A^; is the total number of the lattice sites, and Jo (a;) is the Bessel function of the first kind. 

We next consider H/,. Because there is one link in. Sf and also one in H'^,, the expansion of 
the Baker-Hausdorff formula shows that for the odd terms there appear odd numbers of U, 
combinations. These terms cannot exist in (0| ilj^ |0). We finally get 

(0| 4 |0) =So 5 (-1)^ (0| xHx)xix) |0) 

n=0 (15) 

= ^(-f)*2 E ^.i-mv2er^' 

n=0 

= -if Ji(2y2^, 
with Ji (x) being the Bessel function of the first kind. 

Finally, we consider the last term Hg. Hg can be rewritten as 
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(0| ij; |0) = (0| eM-i0Sf) (^i^ E L\x) j eMiOSf) |0) 

= (0| £ E exp(-i^5/)L(x) ex.p{ieSf) ex.p{-ieSf)L{x) exp(z05/) |0) (16) 

= {0\iE [L{x) + i?(a;) + C{x)] [L{x) + B{x) + C{x)] |0) , 

where 

(0|i?(x)|0)= E i2^(-?)^"^^(-ir-^^E(-i)^ E E mei') 

n=0 ^ ^ 2/ m=0^'=^_2„ (17) 

^+(y)[/2n+l-m([j+)m^(y + 2n - 2m + l)(-l)^"-^"'+"''+^<5j^,^+^' • 



(01 C{x) 10) = E ? ^0 . 2 ^-^^"'^ 

"=1 2/ "»=0m'=m-2n+l (18) 

X+(y)f/2"-"(^7+)'"x(y + 2n- 2m)(-l)2"-2"^+"^'+ij^^^^^,. 
We now only need to calculate B'^{x) and C^(x) in: 

(o| |o) = C E [(o| ^'(^) |o) + (o| C\x) \0) . 

Here we use the fact that L{x) |0) = and (0| B{x)C{x) |0) = 0. The reason for (0| B{x)C{x) |0) = 
is that (0| B{x)C{x) |0) involves (0| x^{y)x{y + 2n-2m + l)x^{y')x{y' + 2n' - 2m) |0). The only 
chance that it will not be zero is that y' = y + 2n — 2m + 1 and y = y' + 2n' — 2m' . This implies 
that 2n — 2m + 1 = 2m' — 2n' . This is not possible because the left hand side is odd and the right 
hand side is even. Therefore, we have 

i E B^x) 10) 4(f )i E ^(V2.)-(-i)"- (:) (::,^) , (19) 



E (01 c\x) 10) = E 2-^j^{^2er{-irOO 

E 2-^^{^ef-{-ir{':) (20) 



2a\ 2 ' Tin) 
n=l 

We wind up with 



- ^(f ) E 72^(V2^)'"(-ir(n) " \ 



(0|ij;|0) 



E 

n=l 



(2n) 



,(^^)2n(_i)„(2n) ri(2n) 



1 /2«-2\ 

2 V n-1 J 



= i(f ){|Jo'(2V2e) - \M2V2e) + I -llo'^'Uodydx [Jf(y) - Ji{y)]}. 



(21) 
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The dimensionless ground state energy now reads 



en = ^= -^moaJo(2v^)-^Ji(2v^) + ^{|j2(2V2^)-iJo(2V20) + | ^^^^ 

-Uo'^''lodydx[jny)-Ji{y)]. 

By = 0, we can determine the parameter 9 for any coupUng constant and bare fermion 



mass. 



4 Results and Discussion 

Based on the unitary transformation and the variational method we addressed in the last sec- 
tion, we have studied the vacuum structure and chiral symmetry breaking in the Schwinger model 
with Kogut-Susskind fermions. The reason we prefer to use the Kogut-Susskind formalism is its 
simplicity. 

The fermion condensate can be calculated by 

^) = {n\ J2 {-irx^{x)x{x) \n) /Ni = -i Jo(2V2^). 

X 

Figure 1 plots (^I/ ^)/(? as a function of in the chiral limit mo = 0. It shows a very good 
scaling behavior although the result we calculated, (^f ^) = —0.252 ± 0.034, is higher than the 
exactly calculated result: 

^)/e = - exp(7)/27r=^/2 ^ _o.i6. 

Our result would be better if we add the interaction with next to nearest-neighbor term to Sf. 
Figure 2 shows (vj/ vj/^ against the bare fermion mass for l/g^ = 1.0, 2.0, and 3.0, respectively. 

For quite a long time, LGT cannot escape the restrict of strong coupling regime. With the 
application of a unitary transformation and the variational method, the awkward situation of LGT 
has been overcome. In a future work, we will use the technique to study the chiral symmetry 
properties and mass spectrum of a two dimensional SU{2) gauge theory coupled to an unfiavored 
Susskind fermion. 

This work was supported by the National Science Council of Republic of China under Grant 
No. NSC 83-0208-M-006-015. We thank Professor Su-Long Nyeo for fruitful discussions. 
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Figure Captions 



Fig.l — ^)/g was plotted as a function of Ijg^ in the chiral limit mo = 0. ranges from 0.5 
to 10, where the symbols square and circle represent our results and the exactly calculated value, 
respectively. 

Fig.2 —{^ was plottted as a function of bare fermion mass rriQa for l/g^ = 1.0,2.0, and 3.0, 
where the symbols cross, square, and triangle represent l/g^ = 1.0,2.0,3.0, respectively. 
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